UMTG-243 



Finite size effects in the XXZ and sine-Gordon models 

with two boundaries 



Changrim Ahn ^ and Rafael I. Nepomechie 



Abstract 

We compute the boundary energy and the Casimir energy for both the spin-^ XXZ 
quantum spin chain and (by means of the hght-cone lattice construction) the massive 
sine-Gordon model with both left and right boundaries. We also derive a nonlinear 
integral equation for the ground state of the sine-Gordon model on a finite interval. 
These results, which are based on a recently-proposed Bethe Ansatz solution, are for 
general values of the bulk coupling constant, and for both diagonal and nondiagonal 
boundary interactions. However, the boundary parameters are restricted to obey one 
complex (two real) constraints. 
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1 Introduction 



The spin- 1 XXZ quantum spin chain and the sine- Gordon quantum field theory on a finite 
interval (i.e., with both left and right boundaries) have applications in statistical mechanics, 
condensed matter physics and string theory, and have therefore been studied intensively, e.g., 
P - ^H] • Much of this work has been restricted to either diagonal boundary interactions - 
[Z] - [ini or to special values of the bulk coupling constant [131 Ellin], because a solution 
of the XXZ chain with general (both diagonal and nondiagonal) boundary terms [5] has 
not been available. A solution of the latter problem for values of the boundary parameters 
obeying a linear constraint has recently been proposed jTHl El E] and confirmed numerically 

m- 

We exploit here this new solution to compute finite-size corrections to the ground-state 
energy of both the XXZ chain and (by means of the light-cone lattice approach [2111121122]) 
the massive sine-Gordon model, in a range of parameter space heretofore not possible. In 
particular, we compute the boundary energy and Casimir energy, and we derive a Kliimper 
- Pearce - Destri - de Vega [^ |21] nonlinear integral equation for the ground state of the 
sine-Gordon model on a finite interval. 

The outline of this article is as follows. In Section |21 we consider the open XXZ quantum 
spin chain with spins. After a brief review of the Bethe Ansatz solution p!Hl IT71 ITHj . 
we compute the ground-state energy, in particular the corrections of order 1 and 1/A^, and 
therefore [26l l27j . the central charge. In Section we turn to the sine-Gordon model 
on an interval of length R. We observe that this model contains an additional boundary 
parameter I7+ — 7_| which has not previously been noted. We analyze the light-cone lattice 
[312011211221 version of this model, which is formally quite similar to the open XXZ chain. We 
determine the relation between the lattice and continuum boundary parameters by matching 
the boundary (order 1) energies in the corresponding models. We then formulate a nonlinear 
integral equation [2S1 121] for the ground state, and give a corresponding formula for the 
Casimir (order 1/-R) energy. In the ultraviolet {R ^ 0) limit, the central charge of the sine- 
Gordon model coincides with that of the XXZ spin chain. Our result for the Casimir energy 
at the free-Fermion point coincides with the result from the TBA approach of Caux et al. 
[TK] . Moreover, we compute the Casimir energy numerically over a wide range of bulk and 
boundary parameters, and track the crossover from the ultraviolet to the infrared regions. 
We conclude in Section jH with a brief discussion of our results and some interesting open 
problems. 
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2 The open XXZ quantum spin chain 



We begin by briefly reviewing the recently-proposed fHl ^1 QBj Bethe Ansatz solution of the 
open spin- 1 XXZ quantum spin chain with both diagonal and nondiagonal boundary terms. 
In terms of the parameters introduced in the latter reference, the Hamiltonian is given by 



N-l 



= l{Yl + ^n^n+l + cosh 1] 

n=l 

+ sinh 1] ^ coth q;_ tanh /5_cr^ + csch a_ sech /?_ ( cosh O^a^ + i sinh 0_af) 

— coth a+ tanh + csch a+ sech /5+ ( cosh 6'+cr^ + z sinh 6'+cr^) |, (2.1) 

where .a^.a"" are the usual Pauli matrices, rj is the bulk anisotropy parameter, a±, [3± , 9± 
are boundary parameters, and is the number of spins. The boundary parameters are 
assumed to satisfy the linear constraint 

a_+/3_ + a+ + /3+ = ±(^_-^+)+r/fc, (2.2) 

where k is an even integer if N is odd, and is an odd integer if N is even. In terms of the 
"shifted" Bethe roots {uj} [1^, the energy eigenvalues are given by 

E = sinh^r7> — -— — -— — -\ — sinh r? (coth a_ + tanh /?_ -|- coth a+ + tanh /?+) 

^ smh(Mj - 2) smh(Mj + 2)2'^ ' 

+ i(iV- l)cosh?7, (2.3) 

and the Bethe Ansatz equations are given by 

sinh(-Uj + f ) \ sinh(2uj + r/) sinh(Mj — f aJ) cosh(Mj — f + /3_) 



sinh(Mj - ^)/ sinh(2Mj - rf) sinh('Uj + ^ - aJ) cosh(Mj + f - /5_) 

sinh(Mj — ^ + «+) cosh(£tj — f + /?+) -i^ sinh(Mj — Uk + r]) sinh(Mj + Uk + r]) 
sinh('Uj + f — tt+) cosh('Uj + ^ — /?+) sinh('Uj — Uk — v) sinh('Uj +11^ — 7]) 

j = l,---,M, (2.4) 
where the number M of Bethe roots is given by 

M = ^{N-l + k), (2.5) 



k being the integer appearing in ()2.2|1 . The case of diagonal boundary terms |2ll3j corresponds 
to the limit j3± —>■ ±00, in which case the constraint ()2.2|1 disappears. 
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We restrict our attention here to the "massless" regime (bulk anisotropy parameter t] 
purely imaginary, with < rj < tc); and therefore, to ensure Hermiticity of the Hamil- 
tonian (jTH), we restrict the boundary parameters a± , 9± to be purely imaginary, and P± to 
be purely real. It is convenient to define new bulk and boundary parameters, 

T] = ifi , a± = ifia± , P± = fib± , 9± = ific± , (2.6) 

where fJ- ,a± ,b± ,c± are all real, with < /i < vr. We use the periodicity a± t— > a± + 27rz of the 
Hamiltonian ()2.H) (and in fact, of the full transfer matrix) to restrict a± to the fundamental 
domain — tt + f < '^m a± < vr + ^, which implies 

^-z/<a±<^ + z/, (2.7) 

where u = - > 1. 

Considering separately the imaginary and real parts of the constraint equation ()2.2j) . we 
see that the boundary parameters must in fact satisfy a pair of real constraints 

a_ + a+ = ±|c_ — c+l + A;, 

b_ + b+ = 0. (2.8) 

The absolute values can be introduced without loss of generality, since the preceding sign is 
arbitrary. 

The energy eigenvalues E depend on the parameters c± only through the absolute value 
of their difference, |c_ — c+|. Indeed, by performing on the Hamiltonian ()2.H1 a global spin 
rotation about the z axis, the parameters c± are shifted by a constant, i.e., c± ^ c± + const. 
In particular, one can eliminate one of these two parameters (say, c+), which results in 
a shift of the other (c_ i— > c_ — c+). Hence, the energy depends on the difference c_ — 
c+. Furthermore, by performing on the Hamiltonian a time- reversal (complex-conjugation) 
transformation, the parameters c± are negated, i.e., c± i— > —c±. Hence, the energy in fact 
depends on |c_ — c+|. 

Let us consider the energy of the ground state of this model as a function of A^, for large 
A^. The leading (order N) contribution, which does not depend on the boundary interactions, 
is well known |2H1- Our objective here is to compute the boundary (order 1) and Casimir 
(order 1/A^) corrections. 
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2.1 Boundary energy 



Let us streamline the notation by defining the basic quantities ^ 

sinh//(A+f) ZTT cosh//(A + f) 



e„(A) 



9nW 



cosh/i(A-f) ■ 



M 



sinh/i(A-f) ' 
The Bethe Ansatz Eqs. fl2.4|) then take the compact form 

^ a >i2iV+l ^ /X ^ e2a_-l{Xj) e2a^-l(Aj) 

gi+2ib^[^j) gi+2ib+[^j) j^^-^ 

J = 1 , , , 

where we have set Uj = fiXj, and we use the new parameters introduced in ()2.(i|l 



(2.9) 



- n ^^l-^i - ^k) 62 (Aj + Afc) 



(2.10) 



We wish to focus here on the ground state with no holes. Hence, we take even, since 
states with N odd correspond to excited states with an odd number of holes. Moreover, we 
take (see Eq. (Q) 



k 



M 



N 



(2.11) 



According to ^H], for this case the Bethe Ansatz solution correctly yields the energy of the 
ground state, and the shifted Bethe roots corresponding to this state are real. However, we 
have subsequently found through further numerical studies of chains with small values of 
A^ that this statement must be qualified: there are regions in the parameter space ()2.7|) for 
which some of the shifted Bethe roots are imaginary (presumably corresponding to boundary 
bound states), or for which the Bethe Ansatz does not yield the ground state. (See Figure 
^) For simplicity, we henceforth restrict the boundary parameters a± to the following four 
regions. 



/ 
// 
/// 
IV 



I < a± < I + z/ 

|<a+<| + z/ & 
1 

2 

i-z/<a+<0 & 



< a_ < 



(2.12) 



I < a_ < i + z/ 



for which our numerical results indicate that the Bethe Ansatz solution does yield the energy 
of the ground state, and the shifted Bethe roots corresponding to this state are all real. 

We remark that (j2.7|) implies that — 2z/ < a_ + a+ — 1 < 2z/; and hence the first constraint 
in Eq. ()2.8|) with k = 1 implies 



|c_ — c+l < 2z/ . 



(2.13) 



^We follow the notations used in [TH] . 
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This condition can always be satisfied, since the Hamiltonian and transfer matrix also have 
the periodicity 6±^ 6± + 2m, which corresponds to c± t— > c± + 2z/. 

In order to compute the energy of the ground state for large N , we first determine the 
density of (real) Bethe roots for this state. To this end, we take the logarithm of the Bethe 
Ansatz Eqs. ()2.10|) and obtain 

MA,)=j, j = l,...,^, (2.14) 
where the counting function h{X) is given by 

KX) = — |(2iV+l)gi(A)+ri(A) + g2a_i(A) — f l+2ib- (A) — f 1+2*+ (A) 



JV 
2 



- Y,M\-\k) + q2{\ + \k)]]. (2.15) 

k=l 

where g„(A) and r„(A) are odd functions defined by 

Q'n(A) = TT + i lne„(A) = 2 tan~^ (cot(nyu/2) tanh(yuA)) , 

r„(A) = ^ln^?„(A). (2.16) 

We have checked numerically that, for the ground state, the right-hand-side of ()2.14|) is 
indeed given by successive integers from 1 to N/2 mi2]- The Bethe roots {A^} can all be 
chosen to be strictly positive. Then, defining A_fc = — A^, we rewrite the last term in ()2.15|) 
more symmetrically as follows: 

iV JV 
2 2 

- fe(A - Afc) + g2(A + Afc)] = -Y, ^2(A - A^) + q^^X) . (2.17) 
fc=i k--^ 

The root density p(A) for the ground state is therefore given by 
1 dh 

P^^^ = NdX 

/>oo 

= 2ai(A) - / dX' a2(A - A') p(A') (2.18) 



+ ^ [ai(A) + 6i(A) + a2(A) + a2a_-i(A) - &i+2i6_(A) + a2a+-i(A) - &i+2ib+(A)] , 

where we have ignored corrections of higher order in 1/A^ when passing from a sum to an 
integral, and we have introduced the notations 

Id sin(n/i) 

«n(A) = 7r-TTQn{X) 



2tt dX TT cosh(2yuA) — cos(n/i) ' 

, ^ d u sinfrau) 

27r dX 71 cosh(2/iA) -|- cos(n/i) 
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The linear integral equation ()2.18p for p(A) is readily solved by Fourier transforms, 



p(A) = 2.(A) + -R{X) 



(2.20) 



where 



-ILUX 



-oo 
oo 



duj e 



2cosh(w/2) 2cosh(7rA) 
^iu^x f sinh((z/ - 2)^/4) cosh(z/^/4) 



[2. 



sinh((z/- 2)cj/2) 



+ sgn(2a+ - 1) 



sinh((z/ - l)cu/2) cosh(cj/2) 
sinh((z/ - |2a+ - 1|)cj/2) 



2sinh((i/ - l)cj/2) cosh(cu/2) 
sinh((l + 2z6+)cj/2) 



2sinh((z/- l)cj/2)cosh(c<;/2) 2sinh((z/- l)cj/2) cosh(cj/2) 
which we have obtained using the results ^ 

sinh {{u — |n|)co'/2) 



sgn[n 



sinh {vuj/2) 
sinh {tujJ /2) 



sinh {vuj/2) 

where = - > 1, and the sign function sgn(n) is defined as 



< |n| < 2z/, 
< 5Re n < z/ . 



(2.22) 
(2.23) 



sgn^nj 



|n| 





n ^ 
n = 



(2.24) 



We have also made use of the fact \2a± — 1| < 2z/, which follows from ()2.7p . 

Having determined the root density for the ground state up to order we now proceed 
to compute the corresponding energy. Recalling the result ()2.3p for the energy in terms of 
the Bethe roots, we obtain 



E 



2tt sin fi 

vr sin /i 
/i 



N_ 

2 



^ai(A,) + ... 



TT sm /i 



2 



^ ai(Aj) -ai(0) 



7 = --^ 
■J 2 



AT 



(iA ai(A) p(A) — ai(0) 



+ . . . 



(2.25) 



where again we ignore corrections that are higher order in and the ellipsis (...) denotes 
the terms in ()2.3p which do not depend on the Bethe roots. Substituting the result ()2.20|) 
for the root density, we arrive at the final result for the ground-state energy 



^Our conventions are 



fWdX, /(A) 



(2.26) 



1 

2^ 



-iuX 



/(w) duj . 
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where 



Ebuik = / rfA ai(A) s(A) + -TV cos/i 

J-oo 2 



= -A^sinV / d\ — — ^ — — - + -A^ cos /z , (2.27) 

J-oo (cosn(2yuA) — cos/i) cosh(7rA) 2 

which is the well-known j2H] result for the bulk (order A^) ground-state energy of the XXZ 
chain; and the boundary (order 1) energy is given by 



71 sin /"oo 

Eboundary = / dX ai(A) [R{X) - (5(A)] 

J-oo 



+ - sin /i (cot /ia_ + i tanhyu6_ -|- cot yua+ + z tanh/i6+) — - COS /i , (2.28) 

where -R(A) is given by 1)2.211) . It should be understood that the boundary parameters obey 
the constraints ()2.8|) with k = 1. In the limit of diagonal boundary terms b± — > ±oo, this 
result for the boundary energy agrees with that of |3j. 

The result ()2.28|) for the boundary energy, which is the sum of contributions from both 
boundaries, implies that the contribution of each boundary is given by 

, _ sin/i /•°° 1 j-sinh((i/-2)cu/4) 1 



boundary ^fi J _^ 2 COsh(cj/2) I 2 sinh(z/Cj/4) 2 

, , sinh((z/ — |2a-|- — l|)u;/2) sinh(ci;/2) cos(6±ti;) 
+ sgn(2a± - 1) — —77 \ , , \ / , , + 



2 sinh((i/ - l)uj/2) cosh(u;/2) sinh(i/u;/2) 

11 , , 

-|- - sin /i cot /la-t — - cos /i . (2.29) 

Indeed, as already noted, the total energy depends on c± only through the combination 
|c_ — c+|. Hence, the left and right boundary energies must be independent of c±. ^ Let us 
now consider the left-right symmetric boundary case, with 0+ = a_ and 6+ = — 6_ (and c± 
are arbitrary, so that a+ = a_ is arbitrary). For this case, we expect that the energy contri- 
butions of the left and right boundaries are equal, E^^^^^y = E^^^^^^^y. Dividing ()2.28p in 
half, we obtain the result ()2.29p . We now argue that this result holds for the general (nonsym- 
metric) case. First, the form of the Hamiltonian 1)2.11) implies that the functional dependence 
of the right boundary energy on the right boundary parameters should be the same as the 
functional dependence of the left boundary energy on the left boundary parameters; i.e., 
boundary = fja+^K) and E-^^^^^y = f{a. , 6„) with the same function /. The boundary 

•^For example, consider the right boundary energy Eboundary- does depend on c+, then it must also 
depend on c_, since the dependence must be of the form |c_ — c+|. But the right boundary energy can 
depend only on the right boundary parameters. Hence, it cannot depend on the left boundary parameter 
c_; and therefore, it cannot depend on c+. 
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energy expressions ()2.29|) evidently satisfy this property. Finally, the left and right bound- 
ary energies must be independent of each other. Hence, having computed E^oundaryi^+ '^+) 
for the left-right symmetric boundary case for arbitrary a+ and it cannot change if we 
change a_ and/or (Although we do not know the Bethe Ansatz when ^ — one 
could in principle do the computations numerically.) Thus, the expression E^oundaryi^+ ' ^+) 
must be correct even for the left-right nonsymmetric case. Similarly, E^oundary^^- '^-) '^ust 
be correct even for the left-right nonsymmetric case. 



2.2 Casimir energy 



The computation of the Casimir (order energy requires considerably more effort. A 

systematic approach based on the Euler-Maclaurin formula and Wiener-Hopf integral 
equations |2H] was developed in [201 for the periodic XXZ chain, and was extended to the 
open XXZ chain with diagonal boundary terms in jS] . Fortunately, the analysis of our system 
of Bethe Ansatz equations ()2.3|) . ()2.4|) is very similar to the one presented in [3]. Hence, we 
shall indicate only the significant differences which occur with respect to this reference, which 
we now denote by I. Using the fact that g„(A) — > sgn(n)7r — /in for A ^ oo, we find that the 
"sum rule" (13.30) becomes 



r°° 1 



2(1+- 



1 , 

— (a_ 

V 



a+ — — ib^ 



where s± = sgn(a-|- — |); and the quantity a (13.34) becomes 



a 



G+(0) 



z6_ 



ib^ - 1] 



(2.30) 



(2.31) 



where, in our conventions, G'+(0)^ = 2(z/ — The Casimir energy (i.e., the contribution 

of order to the ground-state energy) is given by 



Casimir 



CTT sm /i 

24/iiV 

where j25^ I77j the central charge c is given by (13.38) ^ 



(2.32) 



1 - 12a' 



1 - 



6 



V 



[s- + S I 



[a_ + — ib- — ib^ — 1] 



(2.33) 



*In the diagonal limit, the corresponding result is 



c= 1 



6 



- (s_ + s+) - (a_ + a+) 



In particular, the central charge equals 1 for the case a+ + a_ =0 where the boundary fields are real and 
opposite as well as for the case a± = vjl of vanishing boundary fields. In [31, it is implicitly assumed 
that a± > 1/2, in which case s± = 1. 
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Imposing the constraints ()2.8|) with k = 1 gives the final resuh 



c= 1 - 



6 



- (s_ + s+) T \c- - c^ 



(2.34) 



Since the root density should be nonnegative, it follows from the sum rule ()2.3U|) that 
the boundary parameters should obey 



i (1 + s_ + s+) - (a_ + a+ - ib^ - ib+ - I] 



i (1 + s_ + s+) T^\c-- c+l > ,(2.35) 



i.e., ±|c_ — c+l < u{l + S- + s+)/2, which is a further restriction of the constraint ()2.13|) . 

3 The sine-Gordon model with two boundaries 

We turn now to the sine-Gordon quantum field theory on the finite "spatial" interval x G 
\X— , 1, with Euclidean action 



dy / dx A{(p , d^ip) + 



A 



where the bulk action is given by 



dy 



dip, 
dy' 



dip, 
dy' 



x=x+ 



, (3.1) 



A{ip , d^(p) = -{d^ipf + ixbuik cos(/3v?) 



and the boundary action is given by 



i?.(^,^)=/i.cos(-(c,-y,o))±-^^ 



(3.2) 



(3.3) 



This action is similar to the one considered by Ghoshal and Zamolodchikov ^6], except 
that now there are two boundaries instead of one, and the boundary action ()3.3|) contains 
an additional term depending on the "time" derivative of the field. In the one-boundary 
case, such a term can be eliminated by adding to the bulk action ()3.2j) a term proportional 
to dxdyif, which has no effect on the classical equations of motion. However, in the two- 
boundary case, one can eliminate in this way only one of the two 7-1- parameters (say, 7+), 
which results in a shift of the other (7. 1-^ 7^—7+). Notice that this discussion is completely 
parallel to the one for the 9± parameters of the open XXZ chain ()2.1|) . Indeed, we shall argue 
below that these two sets of parameters are related (j3.21|) . 

Let us consider the energy of the ground state of this model as a function of the interval 
length = x+ — x_, for large R. The leading (order R) contribution, which does not depend 
on the boundary interactions, is well known The boundary (order 1) correction is also 
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known fH^]. Our main objective here is to compute the Casimir (order 1/-R) correction, 
and to derive a nonhnear integral equation [23 1211 12^1 for the ground state. We proceed, 
following the analysis of the case of Dirichlet boundary conditions, by considering the 
light-cone lattice [2011^122] version of this model, defined on a lattice with spacing A. This 
lattice model is formally quite similar to the XXZ chain considered in the previous Section, 
the main difference being the introduction of an alternating inhomogeneity parameter ±A. 
The continuum limit consists of taking ^ oo, A ^ 0, and A —>■ oo, such that the length 
R and the soliton mass m (whose relation to fituik is known [32]) are given by 

R = NA, m = ^e'^^ , (3.4) 

respectively. 

In this approach it is evidently necessary to know the relation between all the parameters 
of the lattice model and those of the continuum quantum field theory (j3.H) - (j3.3p . The 
relation between the lattice and continuum bulk coupling constants is well known (see, e.g., 
[ZmU]): = 8(7r - /x) = 87r(z/ - l)/z/, and therefore ^ 

A.^-:^^. ,3.) 

However, corresponding relations for the boundary parameters have been known only for the 
special case of Dirichlet boundary conditions |2] . 

We determine the general relation between the lattice and continuum boundary parame- 
ters in Section im bv matching the boundary energies in the lattice and continuum models. 
We then formulate a nonlinear integral equation for the ground state, and give a corre- 
sponding formula for the Casimir energy. We examine the ultraviolet (i? — 0) limit, and 
also compare our result at the free-Fermion point (A = 1) with that of the TEA approach 
[T3] . Moreover, we compute the Casimir energy numerically over a wide range of bulk and 
boundary parameters, and track the crossover from the ultraviolet to the infrared regions. 

^It should be clear from the context whether the symbol A refers to the value H3.5|l of the bulk coupling 
constant or to the rapidity variable, as in Also, we note that in the bulk coupling constant /j, 

is restricted to the range < fj, < ^, and the Hamiltonian has a coefhcient e, so that the "repulsive" and 
"attractive" regimes correspond to e = +1 and e = — 1, respectively. Here we instead let ^ have an increased 
range < < tt, and consider a single sign of the Hamiltonian, corresponding to the repulsive regime in 
|l()j . Thus, here the repulsive and attractive regimes correspond to the ranges < /i < ^ and ^ < fi < it, 
respectively. In terms oi — 7r//i, these ranges are v > 2 and 1 < ly < 2, respectively. 
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3.1 Boundary energy and boundary parameters 

For the light-cone lattice model, the Bethe Ansatz equations can again be written in the 
logarithmic form ()2.14|1 . except that the counting function is now given by 

hiX) = ^{iV[gi(A + A) + gi(A-A)]+gi(A)+ri(A) + g2a_-i(A)-ri+2ib_(A) 

JV 
2 

+ g2a+-i(A)-ri+2ife+(A)-^[g2(A-Afc) + g2(A + Afc)]}, (3.6) 

fc=i 

which depends on the inhomogeneity parameter A. 

The computation of the ground-state root density to order proceeds as in Section 
12.11 and we obtain 

p(A) = .(A + A) + s{X - A) + li?(A) , (3.7) 
where s(A) and -R(A) are given by (j2.2H) . Moreover, following [HI 122]) the energy is given by^ 

N_ 

E = -i^[ai(A,- + A) + ai(A,--A)] 

i=i 

= -^|Ary""rfAai(A-A)p(A)-ai(A)| , (3.8) 
where A is the lattice spacing. Substituting the result (j3.7j) for the root density, we obtain 



where 



Ebuik = --r rfAai(A-A)[s(A + A) + s(A-A)] 

^ J — CXD 



2ttA sinh(i/co'/2) cosh(co'/2) 

and 

1 



oo 



Eboundary = / ^A tti (A - A) [R{X) - 5(A)] 



oo 



1 f°° cos(tuA) sinh((z/- l)a;/2) 
duj 



2ttA y„^ sinh(z/cj/2) 



R{lj) - 1 . (3.11) 



^We consider explicitly here only contributions to the energy which depend on the Bethe roots. 
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Taking the continuum limit A^— >oo,A— >0,A^oo, keeping the length R and the soliton 
mass m fixed according to ()3.4j) . we obtain (closing the integrals in the upper half plane and 
keeping only the contribution from the pole at = in) 

1 



bulk 



-m^i? cot (z/7r/2) 



and 

Eboundary 



m 

' 2" 



cot(z/7r/4) 



^ ^ cos((i/ — 2s+a+)7r/2) ^ cosh(7r6+) ^ 



(3.12) 



(3.13) 



sin(z/7r/2) sin(z/7r/2) 

where s± = sgn(a-|- — |). 

The same result (j3T2|l for the bulk energy was obtained by a TBA analysis in [Oj. Using 
the relation (|3.5|) between the lattice and continuum bulk coupling constants, one arrives at 
the well-known result |S2] for the bulk energy of the continuum sine-Gordon model. 

The result fl3.13p for the boundary energy, which is the sum of contributions from both 
boundaries, implies (see the corresponding discussion for the XXZ chain at the end of Section 
12. 1|) that the contribution of each boundary is given by 



boundary 



m 

' 2" 



1 , , , 1 cos((z/ — 2s±a±)7r/2) cosh(7r6-i-) 
-- cot z/7r/4) \ , , J H , , , 

2 ^ ' ' 2 sin(z/7r/2) sin(z/7r/2) 



(3.14) 



Comparing with Al. Zamolodchikov's result fH] IT^ for the energy of the continuum sine 
Gordon model with a single boundary 

m i , , , ^ X 1 X X 1 

2 



icos(7r/(2A)) + isin(7r/(2A)) 



cos(r7/A) + cosh('i9/A) 



2cos(7r/(2A)) 

and using again the bulk relation ()3.5|) . we conclude that the boundary parameters of the 
continuum model {ri± , and of the lattice model (a± , b±) are related as follows: ^ 



U5) 



;i + A - 2Aa±) 



(A + l)(- + za±) 



(3.16) 



v± -- 

= Avrfei = (A + l)/3± , 

where we have also made use of ()2.6p . 

Note that the continuum boundary parameters r]± , ■})± in Al. Zamolodchikov's result 
(I3.15|l are those which appear in the Ghoshal-Zamolodchikov boundary S matrix [Hj. Their 
relation to the parameters fx± , ip^ in the boundary action ()3.3|1 is given by [HI [T2j 



cos ( ^{r]± + i^±] 
cos ( f-(r7± - ii^±] 

OTT 



(3.17) 



^For simplicity, we assume in the remainder of this section that a± > 1/2, and therefore, s± = 1. 
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where 



\ sin (g) 



It follows from ()3.16p that the relation between the boundary parameters of the lattice model 
(a± , /5±) and the boundary parameters in the continuum action , ip"^) is given by 

sinh(a± + ^±) = i^te-^^^o ^ 

sinh(a±-/3±) = (3.19) 

For later convenience, we remark that for the left-right symmetric boundary case with a+ = 
a_ and 6+ = — 6_, these relations imply 

^0 = -Vo = ^l2a+{b+) , /i+ = /i- =/ic|sinh/x(6+ + m+)|, (3.20) 

where the function qn{X) is defined in ()2.16|) . 

We still have not discussed the relation between the lattice parameters 6± and the contin- 
uum parameters 7±. We conjecture that these boundary parameters are related as follows: 

7± = -Avrci = i(A + 1)0± . (3.21) 

We perform a check on this conjecture at the free-Fermion point (A = 1) in Section The 
constraints on the lattice parameters ()2.8j) with k = 1 then imply corresponding constraints 
on the continuum parameters 

V- + V+ = T\1- - l+\ + , 

+ = 0. (3.22) 

Finally, let us verify that the first relation in ()3.16|) is correct in the Dirichlet limit. 
Indeed, in terms of the Ghoshal-Zamolodchikov boundary parameters {C,±,k±), which are 
related to the parameters (?7± ,'d±) by 



1 . o 

„2 ' 



COS cosh = —-^ COS ^± , cos^ r^^ + cosh^ -(9^ = 1 + , (3.23) 



the Dirichlet limit corresponds to k± 0, which implies ^± = rj±. On the other hand, for 
this case, the following relation between lattice and continuum parameters is known [7| llUj: 
C± = Cu'^i' ) f ~ (-'- + ^ ~ 2Aa±)|. This result is evidently consistent with ()3.16p . 
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3.2 Nonlinear integral equation and Casimir energy 



We consider now the computation of the Casimir (order 1/-R) energy. Rather than follow the 
Euler-Maclaurin/ Wiener- Hopf approach of Section 1221 we use instead an approach [23 El! 
based on a nonlinear integral equation for the ground-state counting function ()3.6p . which 
is of interest in its own right. 

The derivation of the nonlinear integral equation for the case at hand is similar to the 
case of Dirichlet boundary conditions treated in p. Indeed, following the usual steps, we 
obtain 

A poo j\l J 

— h{X) = N[si\ + A) + s{\-A)] + R{\)- -^^(A-A' + ^e)— ln(l-e-2"*^(")) 

—G{X - A' - 2e)-^ ln(l - e2-'^(^'+*^)) , (3.24) 
2m dX' 

where e is a small positive quantity, G{uj) = 0,2(00)/ (1 + 02(0;)), and s{X) and -R(A) are given 



by ()2.2Hl . Moreover, the energy is given by 

1 r°° dX 

E = E,^ik + E,,^^,,,y--j _s'(A-A + ze)ln(l-e-2-^("-)) 

+ A/ ^s'(A_A-ze)ln(l-e2-''(^+*^)), (3.25) 
A 2m 

where -Eb^ifc and Eboundary are given by ()3.10|) and 1)3.111) . respectively; and a prime on a 
function denotes differentiation with respect to its argument. Integrating ()3.24|) . taking the 
continuum limit as before ()3.4j) . changing to the rescaled rapidity variable 6 = irX, and 
setting f{6) = 2'Kih{9), we finally obtain the nonlinear integral equation 

2? 1"°° 

f{9) = 2imRsmhe + iPbdryiO) + — d9' S^m G{9 - 9' - le) ln(l - e-^(^'+*^)) , (3.26) 

J -00 

where 

^' 27r 2sinh((i/- l)cu/2)cosh(cu/2) ' ^ ^ 

and Pbdry{9) is the odd function satisfying Pldryi^) = 2R{9). Moreover, 

-E P'bulk ~l~ ^boundary ~l~ -^Casimir ; (3.28) 

where Ebuik and Eboundary are now given by ()3.12p and ()3.13|) . respectively; and Ecasimir is 
given by 

m r°° 

Ecas^m^r = '7^ d9 sinh(^ + le) ln(l - e^(^+^^)) . (3.29) 
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3.2.1 Ultraviolet limit 



Let us now consider the ultraviolet limit R ^ 0. Proceeding as in pi, we obtain Ecasimir = 
— C7r/(24i?), where the central charge c is given by 

c=l--^f^Va-7r)2, (3.30) 



9 . 



and 



O- = Pbdry[00) = TT < 1 



U- 1 



2 + s_) — (a_ + 0+ — ib- — ib+ — 1) , (3.31) 



where s± = sgn(a± — |). We conclude that the value of the central charge for the sine- 
Gordon model coincides with the result ()2.33j) . ()2.34j) for the XXZ spin chain. In terms of 
the sine-Gordon parameters 7± ()3.21|) . the central charge is given by 

2 



6 



-(s+ + s_)t 17+-7- 

2 7r 



(3.32) 



3.2.2 Pree-Fermion point 



A dramatic simplification occurs at the free-Fermion point 0^ = Arc, which corresponds ()3.5|) 
to A = 1, or z/ = 2. Indeed, for this value of the bulk coupling constant, the kernel G{6) 
(j3.27p vanishes. It immediately follows from (j3.26|) that f{9) is given by 

f{9) = 2imRsmh 9 + iPbdry{9) . (3.33) 

Let us now rewrite the expression ()3.29|) for the Casimir energy as 

Ecas^mir = -^^ ^ d9 siuh ^ { ln(l - 6^^^+^^) ) - ln( 1 - c^^^^^^^) ) } ; (3.34) 

and then change integration variables 0' = 6* — y + in the first integral, and 9' = 9 + ^ — ie 
in the second integral. Assuming that the resulting contours can then be deformed to the 
real axis, and dropping the primes, we obtain 

/OO 
d9 cosh^ {in (l - e^(^+f )) + In (l - e'/^^"^)) } 

= _£ p d9 cosh^ In (l - e^(^+f )) (l - g-^^^-^)) . (3.35) 

Using (j3.33|) . we obtain 
m /"°° 

Ecasimir = j d9 cosh ^ lu (l + E,{9) e'^-^-^^^ + E2{9) e-4-«™^^^) , (3.36) 

15 



where 

One can show using ()2.21|1 and that e*^''*^'^^^^^ is given (for A = 1) by 



iPbdryid+'-f) 



smh{{9 + iri+)/2) smh{{9 + ir]^) /2) sinh((^ - ^9+)/2) sinh((^ - ^9_)/2) 



(3.38) 



cosh((^ - iv+)/2) cosh((^ - ir]^)/2) cosh((^ + ^9+)/2) cosh((^ + ^9_)/2) 

and e~^^'"^''y^^~^^ is given by the complex conjugate of the above expression. 

This result can now be compared with the result obtained using the TBA approach of 
Caux et al. il5j. One finds that the Casimir energy is again given by ()3.36|) . with (see Eq. 
(58) in dni) 

Ei{9) = tr [K^e) K+(9)) , £2(6) = det [R^O) K+{9)) , (3.39) 
where K±{6) are the crossed-channel boundary S matrices [S] 

K^i9) = r^C- - 0) { ""^'''^ f 1 . (3.40) 

^2 M _sin(e± + 2^) -^e'^±sinh20 ' ^ ^ 



Note that we have included in the boundary S matrices their dependence on the (real) 

dy 



parameters 7±, corresponding to the ^ terms in the boundary action ()3.3p . ^ The scalar 



factors r±{6) are given by 

COS^i 



r±ie) = ^-^(r^i , -ze) a(^^± , -zO) , (3.41) 



where jSBl 



/ \ COS oc . „ „ . 

a(x,u) = . ^ J ^. 3.42 



2cos(f + 1-1) cos(f-f-f) 



Using the relations ()3.23j) to express K±{9) in terms of the boundary parameters r]± ,'d± , 7±, 
we find that the results ()3.37|1 . ()3.39|1 for Ei{9) and E2{0) agree when the boundary param- 
eters satisfy the constraints ()3.22|) . This is a good check on our results ()3.26p . ()3.29p for the 
Casimir energy for general values of the bulk coupling constant, as well as on the conjectured 
relation ()3.2ip between the boundary parameters 9± and 7±. 



^The relation between the 7± parameters in the boundary S matrix H3.40|l and those in the boundary 
action (|3.3|l is not a priori obvious. The fact that these parameters are the same (and, in particular, that 

dip 
dy 



the normalization of the ^ terms in the boundary action is correct) follows from the observation 6 that a 
shift 7+ ^ 7+ + 7 in the boundary 5 matrix implies a corresponding shift B+{ip , ^) '-^ B+{ip j ^) + -^4^ 



in the boundary action. 
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3.2.3 General values of R and u 



For general values of the length R and of the bulk coupling constant z^, the Casimir energy 
cannot be computed analytically. Nevertheless, one can readily solve the nonlinear integral 
equation ()3.26p by iteration and compute the Casimir energy numerically through ()3.29p .^ 
Sample results are summarized in Figures El - El which show the dependence of Cg// = 
— 24:RE Casimir /t^ ou the various parameters. Note that r = mR. In all cases, the computed 
value of Ceff in the ultraviolet region r — >• agrees with the analytical result (|3.3U|) . (|3.3ip . 
Also, as expected, Ce// — > in the infrared region r —>■ oo. Moreover, one can observe the 
crossover in Cg// from the ultraviolet to the infrared regions. 

These graphs are parametrized in part by the boundary parameters a±,b±, in terms of 
which the function R{6) is defined fl2.21|) . Nevertheless, it is straightforward to translate to 
the sine-Gordon boundary parameters using (|2.6p . (|3.19p - (j3.21|) . Indeed, consider Figure El 
for which a+ = a_ and b± = 0; and therefore ()2.8j) . |c_ — c+| = |2a+ — 1|. This corresponds to 

= and /i^. = /i_ = ficSm{7ia^/u); and also |7_ — 7+| = 7r|2a+ — l|/(z/ — 1). Hence, one 
can infer from this graph the dependence of Cg// on yU+ = /i_ or |7„ — 7+|, keeping Lp^ fixed. 
Similarly, for Figure IH a+ = a_ = 1.4 and 6+ = — which implies = —ip^ = ^q2.8{b+) 
and /i+ = /i_ = /ic| sinh fi{b+ + 21.4) |. Hence, one can infer from this graph the dependence 
of Ceff on /i+ = /i_ or ip^ = —fQ, keeping |7_ — 7+| fixed. 

Finally, we remark that the convergence of the iterative procedure which we use to 
numerically solve ()3.26|) depends sensitively on the values of the various parameters. Owing 
to the great number of parameters, we have not attempted to find the entire domain of 
convergence. 

4 Discussion 

We have exploited the recently-proposed [121 HZl UH] Bethe Ansatz solution of the open XXZ 
chain with nondiagonal boundary terms to compute finite size effects in both the XXZ and 
sine-Gordon models, in a range of parameter space previously not possible. Although we have 
focused here exclusively on properties of the ground state, it should be possible, and quite 
interesting, to generalize this work to excited states, with bulk and/or boundary excitations. 
Such a study has recently been made for the case of Dirichlet boundary conditions jSSl- It 
would also be interesting to introduce a "twist" in the nonlinear integral equation to study 
013 perturbed minimal models with boundaries, and to consider applications of our results 

useful trick 34 is to consider H3.26(l with the shift 6 i-^ 9 + ie, and to work in a range of e (typically, 
centered at e 0.3) for which the Casimir energy does not depend on the particular e value. 
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to condensed-matter systems. 

It would be desirable to investigate these models for the full range of boundary parame- 
ters, unhampered by the constraint Indeed, this constraint precludes an investigation 
of the Casimir energy of the sine-Gordon model in the massless scaling limit as a function of 
X = f(v^d ~ V^o)) which is of interest in certain condensed-matter apphcations [T^ ITU [T3]. 
However, finding a Bethe Ansatz solution for this most general case remains a challenging 
open problem. 
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Figure 1: Domain of boundary parameters a±. For the ruled triangular region, the Bethe 
Ansatz does not give the ground state. For the shaded regions (including the special 
case noted in |19J corresponding to the line 0+ + a_ = 1), the Bethe Ansatz does 
give the ground state, but the shifted Bethe roots are not all real. For the blank 
regions (in particular, those labeled I - IV, as in Eq. ()2.12|) ) the Bethe Ansatz gives 
the ground state, and all roots are real. (Based on numerical results for = 4, 
where b± and c± satisfy Eq. ()2.8|) with k = 1. As N increases, the shaded area also 
increases. We conjecture that for A^ 00, only regions I - IV remain unshaded.) 
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c_eff vs. log r 
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Figure 2: Cg// vs. logr, for a+ — — 1.4 , 1.5 , 1.6, with u — 2.2 and 6+ = — 6_ = 1.3 



c_eff vs. a+ 
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Figure 3: Cg// vs. 0+ = a_, for r = 10 ^ , 0.1 , 0.3 , 1.0 , 10, with u = 1.7 and 6+ = 6_ = 
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c eff vs. b+ 




Figure 4: Cg// vs. b+ = for r = IQ-^o ^ iQ-^ , o.Ol , 0.1 , 0.3, with u = 2.7 and 



c_eff vs. log r 
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Figure 5: Ce// vs. logr, for u = 1.9 , 2.0 , 2.1, with a+ = 1.7 , a_ = 1.5 , 6+ = — 6_ 
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